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O-i' Abstract 

We revisit Bs — )• vri^T, -nK* and pK decays within QCDF formalism and examine the 
\ possible effects of a family non-universal Z' boson in these decays. In our evaluations. 



the strong constraints from — B^ mixing on the strength of the left-handed flavor- 
changing d — b — Z' coupling Bjj^j^ ('^ 0.2 x 10"'^) is also included. Numerically, we find 
that a new weak phase 4>^ ~ —50° involved in Sj"^, the negative combination of the 
u — u — Z' and d — d — Z' couplings and/or Z?^^ with larger absolute value are 
crucial to improve the agreement of B{Bs — )■ tt~'K'^) between the SM prediction and the 



^ , experimental measurement. Generally, B{Bs — t- vr K~^) could be reduced by about 

' at most by Z' contributions. Moreover, combining with the recent updated measurement 

on f3s by CDF and DO collaborations, we also find that Z' effects in Bg — )• tt~K~^ decays 
are very important for examining whether the CKM-like hierarchy persists in a family 
non-universal Z' coupling matrix. 
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1 Introduction 



The fruitful running of BABAR, Belle, CDF and DO in the past years provides a very fer- 
tile testing ground for the Standard Model (SM) picture of flavor physics and CP violations. 
Although most of the measurements are in perfect agreement with the SM predictions, there 
still exist some unexplained mismatches, such as the so-called "vr/i', irir puzzles", anomalous 
Bg — Bg mixing phase, large transversal polarization in i? — )■ (pK*, mismatching forward- 
backward asymmetry in i? — K*fi~^fi~ and so on. Recently, some interesting measurements 
of Bg decays have been reported by CDF collaboration at Fermilab Tevatron, which have at- 
tracted many attentions. Especially, CDF collaboration has made the first measurement of 
Bg TT~K+ decay [D |2] 

B(Bg 7r~K+) = (5.0 ± 0.7 ± 0.8) x 10-^ , 

CDF Collaboration (1) 

AipiBg Tv~K+) = (39 ± 15 ± 8)% . 

The recent theoretical predictions for these two quantities based on the QCD factorization 
approach (QCDF) [3], the perturbative QCD approach (pQCD)[^ and the soft-collinear effective 
theory (SCET) [5] read 

B(Bg tv~K+)qcdf = 8.3 x 10"^ , 
^ QCDF Scenario S4 [6] (2) 

Ai;'{Bg ^ 7r~K+)QCDF = 10.9% ; 
BiBg ^ n'K+),QCD = {7.Qtli ± 0.7 ± 0.5) x 10"^ , 

pQCD m (3) 

B{Bg ^ n-K^)scET = (4.9 ± 1.2 ± 1.3 ± 0.3) x 10"^ , ^^^^ ^ 
A%'-{Bg ^ 7i-K+)scET = (20 ± 17 ± 19 ± 5)%. 

One may note that QCDF (Scenario S4) result for B{Bg — )■ 7T^K~^) is significantly larger than 
CDF measurement Eq. ([1]). The pQCD prediction agree with the data in Eq. ([1]) within its large 
theoretical uncertainty, even though its central value also much larger than the measurement 
Eq. ([1]). In SCET, the contributions involving internal charm quark loops are claimed to be 
non-perturbative and only can be determined by fitting to the data. A recent fit analysis of 
charmless B decays using SCET gives the result Eq. (jlj), which agrees with CDF data. All of the 
theoretical predictions for A'^p{Bg — )■ 7r^/i+) agree with CDF data with their large theoretical 
uncertainties. 
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For this possible branching ratio mismatch, many efforts have been done within both SM 
and possible New Physics (NP) scenarios, for example Refs. [9l [101 HH HZj- A recent detailed 
analysis in Ref. [9] indicates that SU{3) and factorization only remain approximately valid if 
the branching ratio for Bg — )■ tt~K^ exceeds its current value by at least 42%, or if a parameter 
^ describing ratios of form factors and decay constants is shifted from its nominal value by 
more than twice its estimated error. Such analysis also implies that, if Bexp{Bs — > tt~ K^) and 
^ persist in the future, large B{Bs — t- -k" K^) mismatch is possibly induced by NP effects. 

A new family non-universal Z' boson could be naturally derived in certain string construc- 
tions [13], models [H] and so on. Searching for such an extra Z' boson is an important mis- 
sion in the experimental programs of Tevatron [15] and LHC [16]. Within a model-independent 
formalism for a family non- universal Z' model presented by Ref. [T7] , fe — )■ s transition involving 
h — s — Z' coupling has been widely studied. It is interesting that such Z' boson behavior is 
helpful to resolve many puzzles in 6 — )■ s transition, such as "vrii' puzzle" [TH |19] and anomalous 
Bs — Bg mixing phase [201 [21]. So, it is worth evaluating the effects of a non-universal Z' boson 
on 5 — J- (i transition, especially the measured B^ — 7r~K~^. Moreover, with the constraints from 
Bg — Bq mixing, B — ttK, hK*, pK and Xsfi~^p~ decays, a CKM-like hierarchy persists in Z' 
coupling matrix (at least between b — s — Z' and b — d — Z' couplings) [21]. Naturally, it is im- 
portant to examine whether such hierarchy is persisted after the constraints from B^ -k^K^ 
are included. 

In addition, the correlated decays Bg — t- iiK* and pK, which also involves b — )■ dqq [q = u, d) 
transition at quark level, will be measured at the running LHC. Such decays may perform an 
important role in testing SM and searching for NP signal. Especially, due to large theoretical 
uncertainty cancellation, the ratio B{Bs — ?■ p^ K^)/B{Bs — ?■ 7i^K~^) is much more suitable for 
testing SM and probing NP [10] . Therefore, we also evaluate Bg — )■ ttK* and pK decays within 
both SM and a family non-universal Z' model. 

Our paper is organized as follows. In Section 2, we revisit Bg — )■ ttK, ttK* and pK decays in 
the SM within the QCDF formalism; we mainly evaluate the effects of annihilation contributions 
which involve endpoint divergency. In Section 3, we perform a fitting on b — d — Z' couplings 
with the constraints from B^ — B^ mixing and B^ tx^K^ decay; the NP effects on the 
other decays are also evaluated. Section 4 is our conclusions. Theoretical input parameters are 
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summarized in the Appendix. 



2 Revisiting Bg ttK, ttK* and pK Decays within QCDF 
Formalism 



In the SM, the effective weak Hamiltonian responsible for b ^ d transitions is given as [2 

10 

V .V*. iT'-.n" -u r'^n"^ -u v,v*. (r-.n? -u r^ns) - v.,v*. ( ' 



rieS — —F= 



h.c, (5) 



where VqbV*^ {q = u,c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements |25] , Cj the Wilson coefficients, and Oi the relevant four-quark operators. 

In recent years, QCDF has been employed extensively to study the B meson non-leptonic 
decays. For example, all of the decay modes considered here have been studied comprehensively 
within the SM in Refs. [6l [11]. The relevant decay amplitudes for Bs — t- ttK, ttK* and pK 
decays within the QCDF formalism are 
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(6) 



where the explicit expressions for the coefficients af = (M1M2) and /3f = /3f (M1M2) can also 
be found in Ref. [6]. The expressions for the Bg — t- 7ri^* and pK amplitudes are obtained by 
setting (nK) {jxK*) and {rcK) — > {pK), respectively. 

It is noted that the QCDF framework contains estimates of some power-suppressed but 
numerically important contributions, such as the annihilation corrections. Unfortunately, in a 
coUinear factorization approach, the calculation of annihilation amplitude always suffers from 
end-point divergence. In the pQCD approach, such divergence is regulated by the parton trans- 
verse momentum kx at the expense of modeling additional k^p dependence of meson distribution 
functions [1], and a large strong phase is found. In Refs. [HI [26], annihilation diagram is studied 
with SCET and parameterized by a complex amplitude. At present, the dynamical origin of 
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these corrections still remains a theoretical challenge. Within the QCDF framework, in a model- 
dependent way, there are two schemes to regulate the endpoint divergence until now. Scheme I: 
In Ref. [6], to probe their possible effects conservatively, the endpoint divergent integrals are 
treated as signs of infrared sensitive contribution and phenomenologically parameterized by 



X 



Jo y ^ 



(7) 



with pA,n < 1 and (^a,h unrestricted. The different scenarios corresponding to different choices 
of pa,h and have been thoroughly discussed in Ref. [6]. Scheme II: As an alternative 
scheme to the first one, one could quote an infrared finite gluon propagator to regulate the 
endpoint divergent integrals, which have been thoroughly studied in Ref. 
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(a) (b) 
Figure 1: The dependence of B{Bs Tr~K~^) on (I)a,h with pa,h = 1 (colored real line), 
2 (colored dashed line). The black horizontal lines are the experimental data with the solid line 
being the central value and the dashed ones the error-bars (Icr). 



In scheme I, with the central value of the input parameters listed in the Appendix, Fig. [T] 
shows the dependence of B{Bs — t- n^K^) on (j)A{H) with pa{h) = 1, 2 and Ph{a) = 0. As 
Fig. W^a) shows, assuming (pA 7^ (pH, the minimal B{Bs — ?■ 7r^A'+) appears at (pn — 0° and 
(j)A — —110°; Assuming 0^ = (p^, it appears at ^ji^ = 0^ ~ —100°, as shown in Fig. [T] (b). 
With the above (t>H,A values, pa,h = 1 and the central values of the other inputs, corresponding 
above two assumptions for (1>a,h, we get B{Bs — t- n^K^) = 9.0 x 10~^, 9.3 x 10^^ respectively, 
which are also significantly larger than CDF data (5.0 ± 1.1) x 10^^. Moreover, corresponding 
to above two assumptions, we find the direct CP asymmetry for — t- 7T~K~^ is 3%, 6% 
respectively, which are also lower than CDF measurement (39 ± 17)% at 2a level. So, following 
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Table 1: The numerical results for B{Bs 7iK,7iK* , pK){xlO-^). 
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Table 2: The numerical results for ^ tt/sT, tt/sT*, p/r)(x IQ-^). 



Exp. SM SM+Z' 
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Scheme II 


SI 


S2 


Ai^iBs - 

Af;{B, - 


+ 7r-K+) 
+ ttO/^O) 


39 ± 17 13.6l^:^ 
- 19.71}°! 
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-16.71^:1 
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+ p-K+) 
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- 9.6t^:l 


18.31^6 
_iq q+16.3 


9.6l^i 
23.71^^ 


9.61^:^ 

1 fi n+3s-^ 



a reasonable choice presented by Cheng and Chua [TT], in our following calculations, we take 
PA,HiPP,PV,VP) = 1 and = -55°(PP), -30°(Pr), -65°{VP). Such choice of the 

values of parameters pa,h and 4>a.h are similar to the so-called favored scenario "scenario S4" 
for Bu,d decays presented by Beneke and Neubert [6], except that some additional moderate 
values of input parameters (such as = SOMeV and 02 = 0.3) are also used in "scenario 
S4" [6]. For estimating theoretical uncertainties, we shall assign an error of ±0.1 to pa,h and 
±5° to 4>A,H- With such Pa,h, 4>a,h values and the other inputs hsted in the Appendix, we 
present the QCDF predictions for Bg — t- ttK, ttK* and pK decays in the third column of 
Table [H [2] and [3 
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Table 3: The numerical results for A^)?(5, tx^Ks) and A'^f{Bs p°A'0)(x IQ-^). 





Exp. 


SM 

Scheme I Scheme II 


SM+Z' 

SI S2 


A^f{Bs - 


> p^K^) — 


_83+i9 -96l?| 

99+35 70+12 


-89tf -88l^^ 
-56+11 -46lt? 



In scheme II, we adopt the infrared finite gluon propagator derived by Cornwall [27], which 
involves a new parameter gluon mass scale nig, to regulate the endpoint divergency. In Ref. [28] , 
we present our suggestion rrig ~ 500 MeV, which is a reasonable choice so that most of the 
observables for B — t- ttK, ttK* and pK decays are in good agreement with the experimental 
data. Furthermore, compared to the available data for — )■ K~^K^ ,Di*^ K , Bg — t- tttt and so 
on, the most recent evaluation [29] also presents that the gluon mass scale is close to 500 MeV. 
So, in this paper, we shall take rrig = 500 MeV. With such a value, the predictions for 
observables of Bg — ttK, hK* and pK decays are listed in the fourth column in Table [H [2] and 
O We find most of the predictions in scheme I and scheme II are consistent with each other 
except for some observables in Bs p^K^, n^K*^ decays. Because the cancellations among a2 
and —al + §03 + ^(^aew ^'i- & ^'^^ ~^ p^K^, 7r°i^*° decays, the effect of annihilation 
contributions is significant. Furthermore, in B ^ PV decay, the annihilation contributions in 
scheme II provide a larger imaginary part than the one in scheme I [28]. So, the two schemes 
present some different predictions for Bg — )■ p^K^, ir^K*^ decays, which will be judged by the 
upcoming LHC-b and proposed super-B experiments. 

Within both schemes I and II, we find that Bsm{Bs ti" K+) is significantly larger than 
the measurement (5.0 ± 1.1) x 10"^. Due to the dominance of the tree contribution ai in the 
amplitude for Bg ti" K+ decay, the contributions from annihilations are tiny in comparison 
to tti. So, the large B{Bs — t- -k" K'^) problem could not be resolved by tuning the annihilation 
contributions. 

A recent detailed analysis [9] in a flavor symmetry framework indicate that the mismatch 
B{Bs — )■ -K^K^) problem persists within the SM if experimental result of B{Bs — t- -k^ K+) would 
not exceeds its current value (5.0 ± 1.1) x 10~^ [H [2] by at least 50% or if the SU(3)-breaking 
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factor defined by 



fKFBOnjmj^) 
fnFB,K{ml) 




(8) 



is not shifted from the value 0.97lo;5i [9l|22] (corresponding to almost exact SU{3), see Refs. [U 
I22j for detail) by more than twice its estimated error. With the central value Fq''^ = 0.24 
obtained by lattice and pQCD calculations, recent QCDF analysis [TT] gets a good result 
B{Bs —7- 7r~A'+) ^ 5.3 x 10^^. However, as mentioned in Ref. [TT], their result ^ = 1.24 
is larger than 0.97 at 3a. Furthermore, F^"^ is predicted by lattice QCD computations at 
small recoil, while unfortunately charmless B decays happen at large recoil. So, if the future 
refined measurement and theoretical study confirm the large B{Bs — ?■ 7r~K~^) problem, it will 
be a significant signal of NP. In our following study, we will examine the effects of a family 
non-universal Z' model for possible solution. 

3 The effects of family non-universal Z' model 

3.1 Brief review of the family non-universal Z' model 

A possible heavy Z' boson is predicted in many extensions of the SM, such as grand unified 
theories, superstring theories, and theories with large extra dimensions. The simplest way to 
extend the SM gauge structure is to include a new U{1) gauge group. A family non-universal 
Z' boson can lead to FCNC processes even at tree level due to the non-diagonal chiral coupling 
matrix. The formalism of the model has been detailed in Ref. [17] . The relevant studies in the 
context of B physics have also been extensively performed in Refs. [TSl [201 EDI 1^ - 

With the assumption of flavor-diagonal right-handed couplings , the Z' part of the effective 
Hamiltonian for b — )■ dqq {q = u, d) transitions can be written as [191 120] 



where Ol{i = 3, 5, 7, 9) are the effective operators in the SM, and ACi the modifications to the 
corresponding SM Wilson coefficients caused by Z' boson. In terms of the model parameters 
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Table 4: The inputs of the Z' couphngs |21] . 



Solutions 


pL 


pi? 

^ud 


^ud 


^ud 


SI 


2.1 ±3.3 


-0.3 ± 1.3 


-0.5 ±0.4 


0.03 ±0.10 


S2 


0.7 ± 1.9 


-0.1 ±0.8 


-0.3 ±0.2 


0.01 ±0.05 



at the Mw scale, ACj are expressed as 



2 ^ ^ 

4 



AC^9,7 = -7TT7^5.^.(5„';^-5,Y), (10) 



where Bh^^ is the effective chiral Z' coupling matrix element. 



q'q 



Generally, the diagonal elements of the effective coupling matrix B^^^ are real as a result 
of the hermiticity of the effective Hamiltonian. However, the off-diagonal one i?J^ may contain 
a new weak phase 0^. Then, conveniently we can represent ACj as 



AC- 



3,5 



2 






VuVt, 


4 









1 



(11) 



where the real NP parameters (^^^'^^^ ^ll,lr ^^^^ defined, respectively, as 



1 , B 



L 



dd ) 1 



/-LL,LR _ ^ I ^db I zjL,R pL,R _ / TyL,R , c) nL,R 

^ ~ Q T/*T/ \ ud ' ^ud — y^uu -T 

cLL,LR _ 1 I ^db I riL,R D^'^ — ( pL,R _ pL,R\ 

S — Q T/*T/ \^ud 1 ^ud — \^uu ^dd J J 

6 VtdVtb 



0^ = kvg[Bi]. (12) 

Firstly, we shall specify the values of the Z' input parameters in our calculations. In 6 — )■ sqq 
transition, the combinations of 6 — s — Z', u — u — Z' and d — d — Z' couplings (^. e. Cf^'^^ 
^LL,LR^ which could be derived through replacing in Eq. (1121) by 5^) have been well bounded 
by the constraints from B — t- txK^*^ and pK decays [11]. After considering the constraints from 
Bs — Bs mixing, which performs constraint on h — s — Z' coupling solely, one may easily extract 
the results u — u — Z' and d — d — Z' couplings -P^j^ and D^f^ [21], which are recapitulated 
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in Table |1U. Such results are the same as the ones in 6 — )■ dqq transition and will be adopted in 
the following analysis. 

As for b — d — Z' coupling, in Ref. [21], we have found the strength of B^f^ is strongly 
suppressed by — B^ mixing. So, in our following fitting, the constraints from i?^ — B^ 
mixing are also considered. The basic theoretical framework for Bq — Bg mixing within a family 
non-universal Z' model have been given by Refs. [201 EI] in detail. In this paper, we take the 
same framework and conventions as Ref. [21]. Moreover, we also take the recent updated fitting 
results of two mixing parameters Cb^ and (pB^ by UTfit collaboration [23], which are defined as 



- {B,\n!}j\B,) Al^^e^'^l'' ' ^ ^ 



and 



Cp, = 0.95 ±0.14, = -3.1° ± 1.7° (14) 



at 68% probability[23]. 



3.2 Numerical analyses and discussions 

With these theoretical formulae and the theoretical input parameters summarized in Appendix, 
we now present our numerical analyses and discussions. Our fitting is performed with that the 
experimental data is allowed within their 1.68a (~ 90% CL.) error bars, while the theoretical 
uncertainties are also considered by varying the input parameters within their respective regions 
specified in the Appendix. In addition, we quote the scheme I to regulate the appearing end- 
point divergences. Since the Z' parameters -Pjfj"^ and have been severely constrained by 
measured Bg — Bs mixing and B — > -kK decays [21] , we demand that they cannot exceed their 
respective ranges listed in Table HJ where SI and S2 correspond to UTfit collaboration's two 
fitting results for Cb^ and (pBs 

With I-BJ'^I = 0.15 X 10"'^ and the central values of theoretical input parameters, Fig.[2]shows 
the dependence of B{B, K+) on 0^, P^f and D^^. We find that B{Bs n-K+) is 



^In ref. |21) . the errors for the numerical results of the Z' couplings are obtained simply with a Gaussian 
distribution treatment of the outputs caused by the uncertainties of input parameters. In this paper, to evaluate 
the conservative effects of the Z' couplings, we would study the whole possible regions of the Z' couplings, which 
are listed in Table E) 
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Figure 2: The dependence of B{Bs ^ K+) on (/.^, P^f and L>X with l^fj 
and the central values of theoretical input parameters. 
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Table 5: The numerical results for the Z' couplings. 



Solutions 


|i?iJ(xlO-3) 


4>A°) 


pL 


pR 




^ud 


SI 


0.16 ±0.08 


-33 ± 45 


0.45 ± 1.65 


-1.2 ±0.4 


-0.70 ± 0.20 


-0.02 ±0.05 


S2 


0.19 ±0.05 


-50 ± 20 


-0.31 ±0.87 


-0.68 ±0.22 


-0.41 ±0.09 


-0.02 ±0.02 



sensitive to the Z' contributions induced by and -D^^, but dull to and D^^. Interestingly, 
as Fig. |2]^b) and (c) shown, B{Bs — > ^"K^) could be reduced by Z' contributions with negative 
P^ and D^^ at 0^ ~ —50°, which possibly presents a solution for the large B{Bs 7i~K~^) 
problem. 

With B{Bs — 7- n^K^), A'^p{Bs —J- n^K^), Cb^ and (pB^ as constraints, corresponding to 
the two solutions for P^f^ and D^f^, the allowed regions for |Pf(,| and (p^ are shown in Fig. [3] 
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Figure 3: The allowed regions for the parameters \B^j^\ and (/>^ in SI (a) and S2 (b) under the 
constraints from B{Bs — s- Tr~K'^) (blue region), A^p^Bg — ^ Tr~K^) (green region), Cb^ (brown 
region), (pB^ (yellow region) and their combination (orange region) respectively. 



and the numerical results are listed in Table O Interestingly, we find that the new fitting result 
(j)^ = -33° ±45° (-50° ±20°) agrees well with the recent evaluation (pj ~ -48° [SI] through the 
constraints from B — t- tttt. As found in Refs. [201 EI], Bd — Bd mixing performs a very strong 
constraint on B^^^, whose magnitude is suppressed to be < 0.24 x 10"'^ in both SI and S2. 
While, the lower bound 0.08(0.14) x 10-^ in S1(S2) for is dominated by B{B, ti~K+). 
Unfortunately, the constraint from y4^p(i?s — ?■ ti^ K'^) is weak due to its large experimental 
error-bar. As shown in Table HI since -P^^ in S2 are restricted to be smaller than their values 
in SI case [21], the lower bound for \B^^\ in S2 is demanded larger than the one in SI. As shown 
in Fig. El^a) and (b), the allowed ranges for both |-Bf(,| and 0^ in S2 are much smaller than 
the ones in SI. So, the solution S2 is much easier to be excluded by the future more precise 
measurements. Due to the constraint from B{Bs — i- ti^ K'^), the flavor-conserving Z' couplings 
P^d^ and D^'J^ are further restricted. Their numerical results are also listed in Table O 

With the values of Z' couplings listed in Table [5] and the other inputs parameters given 
in the Appendix, including their uncertainties, we present our theoretical predictions for the 
observables in the fifth and sixth column of Table [H [2] and [31 We find B{Bs — -k" K^) is 
reduced by Z' contributions. Due to large theoretical uncertainty, its lower bound 5.8(6.3) 
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xlO ^ agrees with experimental measurement (5.0 ± 1.1) x 10 ^ at Icr error level in S1(S2). As 
a favored solution, within the allowed range for Z' couplings, we take 

15,^,1 = 0.24x10-3, 0^ = -55°, pj^, = -1.6(-0.9), D^, = -0.90 (-0.50) , (15) 

in S1(S2). With these Z' values and the central values of the other parameters as inputs, we 
get 

E{Bs = 7.8 (8.5) x 10"^ . (16) 

Comparing with the central value of the SM prediction 9.5 x 10"^, we find B{Bs — t- tt^K'^) 
could be reduced by about 18% (11%) by Z' contributions at most. Such results means that SI 
is more favored for the large B^^{Bs — t- ti" K^) problem than S2. However, such result B{Bs — ?■ 
'K~K'^) = 7.8 X 10~^ in SI is also 2a larger than the experimental result (5.0 ± 1.1) x 10~^. 

Our evaluations are based on the updated UTfit fitting results [23] for Bg — Bg mixing. 
Due to the constraint from Bg — Bg mixing, a CKM-like hierarchy is demanded in family non- 
universal Z' couplings, I.e., \Bdb/Bsh\ ~ \V^JVtb/V^lVtb\ ~ 0.2 [21]. It is noted that both CDF 
and DO collaborations have updated their measurements of Bg — Bg mixing recently based on 
5.2/6"^ and 6.1/6^^ integrated luminosity, respectively. Their results read 

13s = [0.02,0.52] U [1.08,1.55] 68%C.L. CDF [32] , (17) 
0W = -0.76l[^J^ ± 0.02 68%C.L. DO [33] . (18) 

For comparison, DO result Eq. (fTSi) could also be written as 

Ps = 0.38toil±Om 68%C.L. DO. (19) 

Unfortunately, these results have not been combined together yet by the two collaborations 
and by groups like UTfit [23], CKMfitter [31 [35], or HFAG [36]. It is interesting to note that, 
different from their former combined result [0.27, 0.59] U [0.97, 1.30] [32|, SM expectation for Ps 
~ 0.018 agrees with CDF measurement Eq. (|T7|) at ~ la level. DO result in Eq. (fT9l) . which 
is similar to their former result [0.27,0.59] U [0.97, 1.30] [3Z], is also larger than SM prediction 
for Ps at ~ 2(7 level. So, obviously, CDF and DO measurements have significantly difference 
with each other at small Ps region, where is important for constraining the strength of the NP 
contributions. 
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Figure 4: The allowed regions for the parameters ^^^'^ under the constraints (3s = 

[0.02, 0.52] (C-DF) (green region), AMs (pink region) and their combination (blue region). The 
dashed line corresponds to the assumption l-B^I ~ \^db\ ~ ^-^^ ^ 10"^ as a special scenario. 

If /3s ~ 0.02 is confirmed by the coming experimental measurement, CKM-like hierarchy 
in Z' couphngs matrix possibly could be violated, which implies that \Bsb\ could be as small 
as \B(ih\. To confirm such a guess, we revisit Bg — Bg mixing within the family non-universal 
Z' model following the same track as Ref. [H]. With /3, = [0.02, 0.52] (CDF Eq. ([17])) and 
AMs = 17.77 ± 0.12 [36] as constraints and all of the central values of the parameters listed 
in appendix as input, the allowed regions for the new parameters [-B^l and 0f are shown by 
Fig. m So, from Fig. HJ one may easily find that the relation \Bsb\ ~ j-B^bj ~ 0.2 x 10~^, which 
is smaller than the former result ~ 1 x 10~^ [21] by a factor about 1/5, is allowed by the CDF 
updated measurement on Ps- Due to the constraints from B — )■ irK, pK and ttK*, the values of 
the combining parameters (^^^'^^ and C,s^'^^, whose form could be derived through replacing 
Sfft in Eq. ([H]) by 5^, are unchangeable [HIISI]. So, with \Bsb\ ~ \Bdb\ ~ 0.2 x 10^^ as input, 
after reevaluating the constraints from B — t- ttK et al. (as Refs. [HI |2T] do), one may easily 
find that the results for [-P^j'^l and |-D^^^| could be larger than the ones in Table |4] by a factor 
5. So, clS cl possible special scenario, we take 

[5,^,1 = 0.16 X 10-^ 0^ = -55°, P,^, = -1.2x5, D^, = -0.70 x 5, (20) 

for SI. With such Z' coupling values and the central values of other input parameters, we 
get B{Bs — 7- 'K~K^) = 5.1 X 10^^ which is perfectly consistent with the experimental data 
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(5.0 ± 1.1) X 10"^. Interestingly, as found in Ref. [IT], \Bst,\ ~ \Bdb\ is also needed for resolving 
the "tttt puzzle". It also should be noted that our above analysis is based on the data Ps ~ 0.02 
is allowed, which is also needed to be confirmed by the future refined measurements of Bg — Bs 
mixing system at LHC-b and Tevatron. 

4 Conclusion 

In this paper, we have revisited Bs — ?■ irK, irK* and pK decays within QCDF formalism. We 
find that the large B*^{Bs — t- ti^K'^) problem could not be resolved by tuning the annihilation 
contributions. In order to pursue possible solutions for the large B^^{Bs — )■ ti" K^) problem, 
we evaluate the effects of a family non-universal Z' boson. Bd — Bd mixing is used to constrain 
the h — d — Z' coupling. Our main conclusions are summarized as: 

• The magnitude ioi h — d — Z' coupling is suppressed to ~ 0.2 x 10^^ by Bd — Bd 
mixing. Based on the constrained u — u — Z' and d — d — Z' couplings by Bg — Bg 
mixing and B — ?■ txK decays in Ref. [21], we find -P^^ and D^'^ are further restricted by 
B{Bs ^ TX-K+), especially P^, and D^,. 

• A new weak phase 0^ ~ —50°, negative Pj^ and/or D^^ with larger absolute value are 
helpful to improve the agreement of B{Bs — )■ 7r~A'+) between the SM prediction and the 
experimental measurement. Compared the SM result, B{Bs — )■ -k^ K'^) could be reduced 
by about 18%(11%) at most in SI (S2). SI is more favored by the large tx~ K+) 
problem than S2. 

• The experimental measurement of [ig is very important to fix whether a CKM-like hierar- 
chy is held in a family non-universal Z' couplings matrix. If the lower bound 0.02 for fi^ 
measured by CDF is confirmed, \Bsb\ would be as small as \Bdb\- Within such a scenario, 
B{Bs — 7- 7r^/i+) could be reduced to 5.1 x 10^^ which agrees with the experimental data 
(5.0 ± 1.1) X 10-^. 

The refined measurements for the Bg non-leptonic decay (especially Bg — t- -k" K^) and 
Bg — Bg mixing in the LHC-b will provide a powerful testing ground for the SM and possible 
NP scenarios. Our analysis about the Z' effects on the observables for Bg — )■ ttK, ttK*, pK 
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decays and Bg — Bg mixing are useful for probing or refuting the effects of a family non-universal 
Z' boson. 
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Appendix: Theoretical input parameters 

For the CKM matrix elements, we adopt the UTfit collaboration's fitting results 

p = 0.132 ±0.02 (0.135 ±0.04), = 0.367 ± 0.013 (0.374 ± 0.026), 

A = 0.8095 ±0.0095 (0.804 ±0.01), A = 0.22545 ± 0.00065 (0.22535 ± 0.00065). (21) 

The values given in the bracket are the CKM parameters with assumption of the presence of 
generic New Physics, which are used in our evaluations when the Z' contributions are included. 
As for the quark masses, we take the current quark masses 

= 27.4 ±0.4 [39], mJ2GeV) = 87±6MeV [39], m J fric) = 1 .27^0 u GeV M 
rnb{rfib) = 4.20^2;^^ GeV [IQ] , mt{mt) = 164.8 ± 1.2 GeV [IQ] , (22) 

where mq{fi) = {ffiu + ffidjifi) /2, and the difference between u and d quark is not distinguished. 
The other one is the pole quark mass. In this paper, we take [iQ], HI] 



ITT'u = = TTT's = 0, = l-6lto;i2 GcV, 

-1+0.19 , 



mb = 4.79l^;i| GeV, mt = 172.4 ± 1.22 GeV. (23) 
and the decay constants 



Jb, = (231 ± 15) MeV [12] = (130.4±0.2) MeV [10] = (155.5±0.8) MeV [m 
fx- = (217 ± 5) MeV [13] , {2.2GeV) = (156 ± 10) MeV 
fp = (209 ± 2) MeV [13] , fJ{2.2GeV) = (147 ± 10) MeV 
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As for the B-meson lifetimes, tb^ = 1.470 ps [ID] is used. We take the heavy-to-hght transition 
form factors 

A§'^^*{0) = 0.360 ± 0.034 [13] , F^^^^{0) = O.SOtofs M ■ (24) 



References 

[1] T. Altonen et al. (CDF Collaboration), Phys. Rev. Lett. 103 (2009) 031801, 
larXiv:0812.4"27l] [hep-ex]. 

[2] CDF Collaboration, "Measurement of branching fractions and direct CP asymmetries of 
— )■ h^H'~ decay in lfb~^" and the updated results on April 10, 2008 may be found on 



http://www-cdf.fnal.gov/physics/new/bottom/bottom.html [Note 8579vl] 



[3] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda , Phys. Rev. Lett. 83 (1999) 
1914 |hep-ph/99053l2] ; Nucl. Phys. B591 (2000) 313 | hep-ph/0006T42] . 



[4] Y. Y. Keum, H. N. Li and A. I. Sanda, Phys. Lett. B 504 (2001) 6 |hep-ph/0004"004| ; 



Phys. Rev. D 63 (2001) 054008 |hep-ph/0004173| 



[5] C. W. Bauer, S. Fleming and M. Luke, Phys. Rev. D 63 (2001) 014006 |hep-ph/0005'275] ; 



C. W. Bauer, S. Fleming, D. Pirjol and I. W. Stewart, Phys. Rev. D 63 (2001) 114020 
|hep-ph/0011336] ; C. W. Bauer and I. W. Stewart, Phys. Lett. B 516 (2001) 134 
|hep-ph/010700l1 ]; C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. D 65 (2001) 



054022 [hep-ph/0109045] . 
[6] M. Beneke and M. Neubert, Nucl. Phys. B 675 (2003) 333 |hep-ph/0308039] . 



[7] A. Ali, G. Kramer, Y. Li, C. D. Lu, Y. L. Shen, W. Wang and Y. M. Wang, Phys. Rev. D 
76 (2007) 074018 |hep-ph/0703T62] . 



[8] A. R. WiUiamson, J. Zupan, Phys. Rev. D 74 (2006) 014003, Erratum-ibid. D 74 (2006) 
03901 |hep-ph/06012l4] . 



17 



[9] C. W. Chiang, M. Gronau, J. L. Rosner, Phys. Lett. B 664 (2008) 169, larXiv:0803.3229l 
[hep-ph] . 

[10] G. H. Zhu, JHEP 1005 (2010) 063. larXiv:1002.4518l [hep-ph], 

[11] H. Y. Cheng and C. K. Chua, Phys. Rev. D 80 (2009) 114026. laDGv:0910.5237l [hep-ph]. 

[12] Y. G. Xu, R. M. Wang and Y. D. Yang Phys. Rev. D 79 (2009) 095017, larXiv:0903.0"256l 
[hep-ph]; R. M. Wang and Y. G. Xu, Phys. Rev. D 81 (2010) 055011. 

[13] G. Buchalla, G. Burdman, C. T. Hill and D. Kominis Phys. Rev. D 53 (1996) 
5185 |hep-ph/9510376]; G. Burdman, K. D. Lane and T. Rador, Phys. Lett. B 514 



(2001) 41 |hep-ph/0012073j ; A. Martin and K. Lane, Phys. Rev. D 71 (2005) 015011 
|hep-ph/0404107| . 

[14] E. Nardi, Phys. Rev. D 48 (1993) 1240 |hep-ph/9209"223l ; J. Bernabeu, E. Nardi and 

D. Tommasini, Nucl. Phys. B 409 (1993) 69 |hep-ph/ 93 06251] ; V. D. Barger, M. S. Berger 
and R. J. Phillips, Phys. Rev. D 52 (1995) 1663 |hep-ph/9503"204| ; M. B. Popovic and 

E. H. Simmons, Phys. Rev. D 62 (2000) 035002 |hep-ph/000l"302] ; T. G. Rizzo Phys. Rev. 
D 59 (1999) 015020 [hi^ph/9806397] . 



[15] M. S. Carena, A. Daleo, B. A. Dobrescu and T. M. P. Tait, Phys. Rev. D 70 (2004) 093009 
|hep-ph/0408098] . 

[16] T. G. Rizzo. |hep-ph/0610104| ErXiv: 0808. 19061 [hep-ph]. 

[17] P. Langacker and M. Pliimacher, Phys. Rev. D 62 (2000) 013006 [hep-ph/Q0012Q4] . 



[18] V. Barger, C. W. Chiang, P. Langacker and H. S. Lee, Phys. Lett. B 598 (2004) 218 
|hep-ph/0406126| . 

[19] Q. Chang, X. Q. Li and Y. D. Yang, JHEP 0905 (2009) 056. larXiv:0903.0275l [hep-ph]. 

[20] V. Barger, L. L. Everett, J. Jiang, P. Langacker, T. Liu and C. E. M. Wagner, JHEP 0912 
(2009) 048, larXiv:0906.3745l [hep-ph]; Phys. Rev. D 80 (2009) 055008, larXiv:0902.4507l 
[hep-ph] . 



18 



[21 
[22 

[23 

[24 

[25 

[26 
[27 

[28 
[29 



Q. Chang, X. Q. Li and Y. D. Yang, JHEP 1002 (2010) 082. larXiv:0907.4408l [hep-ph]. 

A. Khodjamirian, T. Mannel and M. Melcher, Phys. Rev. D 68 (2003) 114007 
|hep-ph/0308297] ; Phys. Rev. D70 (2004) 094002 |hep-ph/04Q7226l . 



M. Bona et al, larXiv: 0906. 09531 [hep-ph]; M. Bona et al. (UTfit Collaboration) 



larXiv:0803. 06591 [hep-ph]: online update at: |http://www.utfit.org/UTfit/Results 



G. Buchalla, A. J. Buras, and M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125 
|hep-ph/95 12380] . 

N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; M. Kobayashi and T. Maskawa, Prog. Theor. 
Phys. 49 (1973) 652. 



T. Feldmann and T. Hurth, JHEP 11 (2004) 037 |hep-ph/0408"l88] . 



J. M. Cornwall, Phys. Rev. D 26 (1982) 1453; J. Papavassiliou and J. M. Cornwall, Phys. 
Rev. D 44 (1991) 1285. 

Q. Chang, X. Q. Li and Y. D. Yang, JHEP 0809 (2008) 038. larXiv:0807.4295l [hep-ph]. 
C. M. Zanetti and A. A. Natale. [arXiv:1007.5072l [hep-ph]. 



[30] K. Cheung, C. W. Chiang, N. G. Deshpande and J . Jiang, Phys. Lett. B 652 (2007) 
285 |hep-ph/0604223] ; J. H. Jeon, C. S. Kim, J. Lee and C. Yu Phys. Lett. B 636 (2006) 



270 |hep-ph/0602T56l ; C. H. Chen and H. Hatanaka, Phys. Rev. D 73 (2006) 075003 
|hep-ph/0602"l40] ; R. Mohanta and A. K. Giri, larXiv:0812. 18421 [hep-ph]; J. Hua, C.S. 



Kim and Ying Li, Phys. Lett. B 690 (2010) 508 |hep-ph/0602156] , larXiv:1002.2532 [hep- 
ph]. 

[31] Q. Chang, X. Q. Li and Y. D. Yang. larXiv: 1003.60511 [hep-ph]. 

[32] CDF collaboration, "Updated Measurement of CP Violating Phase /^f^'^'^". Public Note 
10206, May 10, 2010. 

[33] DO collaboration, "Updated Measurement of the CP- Violating Phase (pi^'^'^ Using Flavor- 
tagged Decay 5° ^ J /ixp ", Note 6098-CONF, July 22, 2010. 



19 



[34] CKMfitter Group (J. Charles et al.), Eur. Phys. J. C 41 (2005) 1 |hep-ph/ 0406184] , up- 



dated results and plots available at: http://ckmfitter.in2p3.fr 



[35] A. Lenza et.al. (CKMfitter Group). laFXiv: 1008. 15931 [hep-ph]. 

[36] E. Barberio et al. (Heavy Flavor Averaging Group). larXiv:0808. 12971 [hep-ex], and online 



update at: http://www.slac.stanford.edu/xorg/hfag 



[37] CDF/DO ATsj/Sg Combination Working Group, "Combination of DO and CDF Results on 
AF, and the CP- Violating Phase ^f^'^'^", CDF Public Note 9787, DO Note 5928-CONF, 
July 22, 2009. 



[38] M. Bona et al. (UTfit Collaboration) JHEP 0507 (2005) 028 |hep-ph/0501199) ; JHEP 



0601 (2006) 081 |hep-ph/0606167l ; online update at: |http://www.utfit.org/UTfit/Results 



Q. Mason et al. (HPQCD Collaboration), Phys. Rev. D 73 (2006) 114501 
|hep-ph/0511160l . 

[40] C. Amsler et al. (Particle Data Group), Phys. Lett. B 667, 1 (2008). 

[41] Tevatron Electroweak Working Group, CDF Collaboration and DO Collaboration, 
larXiv:0808.1089l [hep-ex] . 

[42] Elvira Gamiz et al. (HPQCD Collaboration ), Phys. Rev. D 80 (2009) 014503, 
larXiv:0902.18T5l [hep-lat]. 

[43] P. Ball and R. Zwicky, Phys. Rev. D 71 (2005) 014015 |hep-ph/0406232] ; 71 (2005) 014029 
|hep-ph/0412079] ; Phys. Lett. B 633 (2006) 289 |hep-ph/05 10338] . 



[44] G. Duplancic and B. Melic, Phys. Rev. D 78 (2008) 054015. E^v:0805.4170l [hep-ph]. 



20 



